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Abstract 

This paper is a review of recent results on integrable nonholonomic geodesic flows of 
left-invariant metrics and left- and right-invariant constraint distributions on compact 
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1 Introduction 

This paper is a review of recent results on integrable flows on compact Lie groups under 
nonholonomic constraints. We mostly follow papers |24l 1301 IMl I2()j . trying to present their 
results within a unified framework. Furthermore, some new examples of integrable nonholo- 
nomic systems are given. 

1.1 Nonholonomic Geodesic Flows 

We start with basic definitions and settings. Let {Q,ds^) be n-dimensional Riemannian 
manifold Q with a nondegenerate matric ds^and a Levi-Civita connection V, Z? be a non- 
integrable /c-dimensional distribution on the tangent bundle TQ. A smooth path j{t) e 
Q, t £ A is called admissible (or allowed by constraints) if the velocity 7(t) belongs to £'-y(f) 
for all t € A. There are two approaches to define geodesic lines among admissible paths: by 
induced connection as "straightest" lines and by the variation principle as "shortest" lines. 
We shall deal with the first approach which arises from mechanics. 

The admissible path "/{t) is called a nonholonomic geodesic if it satisfies the equations 

T(V^(t)7W) = 0, (1.1) 

where n : TgQ — > Dg, q € Q is the orthogonal projection. 

Equivalently, we can introduce the Lagrangian function / ~ ^{Kq,q), where K is the 
metric on Q also regarded as a mapping K : TQ ^ T*Q. Let q — (qi, . . . , (?„) be some 
local coordinates on Q. The trajectory of the system q{t) that satisfies the constraints is a 
solution to the Lagrange-d'Alambert equations 

One can also write the Lagrange-d'Alambert equations as a first-order system on the 
{n + A;)-dimensional constraint submanifold A4 = K{D) of the cotangent bundle T*Q. Let 
D be locally defined hy p = n — k independent 1-forms a* 

Dg^iCe TgQ, {al, = E = 0, J - 1, . . . , p}. 

i 

Then Ai is locally given by the equations {ag, K~^p) = 0, i = 1, . . . ,p. Let pi = dl/qi, 
I = 1, . . . , n be momenta which together with q provide canonical coordinates on T*Q. Let 
hiq^p) = \{p, Kg^p) be the Hamiltonian function (the usual Legendre transformation of 
L). The equations 1)1.2(1 are equivalent to 

dh{q,p) ,^.jf. . dh{q,p) 
Pi = '^''"~dp — ' i^^,---,n, (1.3) 

where Lagrange multipliers are chosen such that the solutions {q{t),p{t)) belong to Ai. 

As for the Hamiltonian systems, the Hamiltonian function is always the first integral of 
the system. There is also a nonholonomic version of the Nocther theorem (see |38l 1^ 1211 1^ ) . 
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The Noether theorem. Suppose that a Lie group © acts on the configuration space 
Q and that the action is naturally extended to TQ and T*Q. The momentum mappings 
*z : TQ ^ 0* and 4-* : T*Q ^ g* are defined by 

^M.m = = {K,q,^Q), ^*{qMO = iv.iQ). (1.4) 

where is the vector field on Q associated to the action of one-parameter subgroup exp(^^), 
e e = Tjd©. 

Theorem 1.1. Assume that is a section of the distribution D and the one-parameter 
subgroup exp(t^) preserves I (or h). Then is the first integral of the system 1^1. or 

equivalently, is the first integral of 

Invariant measure and integrabiUty. The equations are not Hamiltonian. This 
is why it is still not clear how to define the notion of complete integrability for nonholonomic 
systems (see 0]). However, in some cases they have an invariant measure, a rather strong 
property, which puts the system close to Hamiltonian systems. In particular, if apart from 
the Hamiltonian there exist dim —3 additional independent integrals, then, by the Euler- 
Jacobi theorem, the solutions of (|1.3|) can be found by quadratures. 

The importance of an invariant measure for integrability of nonholonomic systems was 
indicated by Kozlov in jSH], where various examples were discussed (see also (2|). Namely, 
consider a non-Hamiltonian system 

i; = /(a;), x G R", (1.5) 

having an invariant measure ij.{x) dx and m — 2 first integrals Fi (x) , . . . , Fm (x) . If the latter 
are independent on the invariant set Mc = {x E K™, Fi{x) — Ci, i — 1, . . . ,m — 2}, then 
Mc is a two-dimensional submanifold and the fiow on Mc has also an invariant measure. 
Then, accoring to the Euler-Jacobi theorem, solutions of (|1.5|l lying on Mc can be found 
by quadratures. Moreover, if Lc is a compact connected component of Mc and f{x) ^ on 
Lc-, then, if orientable, Lc is diffeomorphic to two-dimensional torus. 

According to Kolmogorov's theorem on reduction of differential equations with a smooth 
invariant measure on a torus (|37|). one can find angular coordinates (pi, ip2 on Lc, in which 
the reduction of equations 1)1.5(1 takes the form similar as in the Liouville theorem: 



where fii, fl2 depend on the constants of motion ci, . . . , Cm-2 only and $ is a smooth positive 
27r-periodic function in ipi,ip2, the density of the induced invariant measure on Lc- 

Therefore, it is natural to call the system 1(1. 5|l completely integrable if it can be inte- 
grated by the Jacobi theorem; or, more generally (see |54l l55j V if the phase space is almost 
everywhere foliated by invariant tori T'^{ipi, . . . , ipk} with the dynamics of the form 

V^l = -TT 7, '^fc = -j-T (1-6) 

<P((pi,...,ipk) <P(ipi,...,ipk) 

The above definition of complete integrability is slightly different from the definition of 
complete integrability of non-Hamiltonian systems given in [SJ I60| . Namely, here we have 
quasi-periodic motions after the time substitution dr — ^~^((p)dt. 

The existence of an invariant measure for smooth dynamical systems and for a class 
of nonholonomic systems with symmetries is studied in |40| and |59| . respectively. Var- 
ious mechanical examples with an invariant measure can be found in |12j . The authors 
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of the paper [541 155| constructed nonholonomic systems on unimodular Lie groups with 
right-invariant nonintegrable constraints and a left-invariant metric (so called LR systems) , 
and showed that they always possess an invariant measure, whose density can be effectively 
calculated. In particular, the motion of a rigid body around a fixed point under a nonholo- 
nomic constraint (projection of the angular velocity to the fixed vector in space is constant) 
is described by an integrable LR system (jMs)- Similar integrable problems on Lie groups 
with left-invariant constraints are studied in |26l 1801 m] . Also, an important example of an 
integrable nonholonomic mechanical system, the problem of rolling of a homogeneous ball 
on a surface of revolution (the Routh problem), was treated in detail in |771IS7| . 

1.2 Chaplygin Systems 

Another approach to the integrability of nonholonomic systems is based on their reduction 
to a Hamiltonian form after an appropriate time rescaling. First, following f 36| and [Hj, let 
us recall some basic facts about the Chaplygin systems. 

Let {Q,l,D) be a nonholonomic system with a Lagrangian I of the natural mechanical 
type, with kinetic energy that correspods to the metric ds^ and the potential function v. 
Assume that there is a bundle structure it : Q N with the base manifold N and let the 
map TT be a submersion, such that TqQ = Dq(BVq for all q. Here Vq is the kernel of T^tt called 
the vertical space at q. Then the distribution D can be seen as a collection of horizontal 
spaces of the Ehresmann connection associated to t: : Q —^ N. Given a vector Xq G TqQ, 
there is a decomposition Xq = X^ + X^, where X^ G Dq, X^ G Vq. The curvature of the 
connection is the vertical valued 2-form B on Q defined by 

BiXq,Yq) = -[X^,YX 

where X and Y are smooth vector fields on Q obtained by extending of Xq and Yq. 

By applying the Ehresmann connection the Lagrange-d'Alambert equations H1.2|l can 
be represented in the form (see 0) 

/die d die \ fdl,. ,\ ^ 

(^^^^,.)^(^,B(.,.)), forall ,eDq, (1.7) 

where lc{q,q) — l{q,q'^) is the constrained Lagrangian. 

Now, suppose that tt: Q^N = Q/& is a. principal bundle with respect to the left action 
of a Lie group 0, and I? is a principal connection, i.e., Z? is a ©-invariant distribution. Let 
the Lagrangian I be also (S-invariant, i.e., 25 acts by isometrics on Riemannian manifold 
{Q,ds'^) and is a 0-invariant function. Then the constrained Lagrangian Ic induces a 
well defined reduced Lagrangian L : TQ — > M via identification TN Ri D/&. The reduced 
Lagrangian L is of the natural mechanical type as well. Its kinetic energy is given by metric 
dsjj and its potential energy will be denoted by V. 

Under the above assumtions, equations (|1.7|l are ©-invariant and induce reduced Lagrange- 
d'Alambert equations on the tangent bundle TN, 

Here S is semi-basic two- form given by the right hand side of (|1.7|) and q = (gi, . . . , g^) 
are some local coordinates on the base space A^. From 1)1.4(1 we see that E depends on the 
curvature of the connection D and on the momentum mapping <i>; . 

The system {Q,l,D,&) is referred to as a (generalized) Chaplygin system (see [361 
as a generalization of classical Chaplygin systems with Abelian symmetries ^7] . 
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Remark 1.1. Note that horizontal and vertical spaces do not need to be orthogonal with 
respect to the metric ds^. In fact, if D is ds^-orthogonal to the leaf of ©-action, then D 
will be an invariant submanifold of the nonconstrained geodesic flow of the metric ds^, and 
the right hand sides of equations H1.7|) will be zero. In this case, dsj^ coincides with the 
submersion metric induced from ds^ . 



Chaplygin's reducing multiplier. Let pi = dL/dqi, i ~ l,...,fc be momenta, gij 
the metric tensor of dsj^ and g^^ the dual metric on T*N. Then the reduced Lagrangian 
has the form L{q,q) — gijCii<ij — V{q). We also introduce the Hamiltonian function 
H{q,p) — \ 'Yl,9^''PiPj + ^{l)- The reduced system (|1.8|l can be rewritten as a first-order 
dynamical system on T*N: 

OH OH 
= — , pi^~- hnj(q,p), i = l,...,fc. (1.9) 

opi dqi 

The functions Hi are quadratic in momenta and can be regarded as non-Hamiltonian per- 
turbations of the equations of motion of a particle on N. 

Let 17 = ^ dpi A dqi be the standard symplectic form on T*N. The equations 1)1. 9|l have 
an invariant measure /O'^ if J2i {^^^ + ^^^^Op ^'^'^ ) ^ ^i^^e the finction / depends 
only on the coordinates q, this is equivalent to condition 

d{\nf) + a^Q, (1.10) 

where the one-form a is given by f^|(j=gp — (q^j i)- 

Remark 1.2. The paper [HO] (see also ^31) contains a nontrivial observation about the 
density of the invariant measure, which in our terms reads as follows. Suppose that system 
(II. 9|) has an invariant measure with density f{q,p) in the case of absence of potential iV{q) = 
0). Then one can check that the function fo{q) = f{q,0) is also a solution of (|1.10() . In 
other words, if the reduced system (|1.9|l has an invariant measure for V" = 0, one can take 
this measure to be of the form f{q)^^. Then, since (I1.1U|I does not depend on the potential, 
the reduced system (11.911 has the same invariant measure in the presence of a potential field 
V{q) as well. 

Now consider time substitution dr = M{q)dt, where M{q) is a differentiable nonvanishing 
function on Q, and denote q' — dq/dr. Then we have the following commutative diagram 

TN{,.j} TN {,.,'} 



'1 1^ 



T*N{q,p} T*N{q,p}. 

The Lagrangian and Hamiltonian functions in the coordinates {q, q'} and {q,p} take the 
form 

«') = ^ E -^'aM, - ViQ), H*{q,p) = i E -^9''P^PJ + V{q). 

There is a remarkable relation between the existence of an invariant measure of the 
reduced system (|1.9|l and its reducibility to a Hamiltonian form (see j26j ) . 

Theorem 1.2. 1). Suppose that after the time substitution dr — J\f{q)dt the equations 
M.9\) become Hamiltonian, 

Op^ Oq^ 
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Then the function f{q) — J\f{q)^ ^ satisfies the equation fLlU^) . i.e., the original 
system ^1.91) has the invariant measure with density f{q). 

2). For k — 2, the above statement can also he inverted: the existence of the invariant 
measure with the density JV{q) implies that in the new time dr = J\f{q)dt, the system 
il.yp gets the Hamiltonian form 

In nonholonomic mechanics the factor J\f is known as the reducing multiplier, item 2) of 
this theorem is referred to as Chaplygin's redueihility theorem (see |16l I17| or section III. 12 
in Notice that for k > 2, the multipher Af{q) and the density of the invariant measure 

of system (|1.9ll do not coincide. Also, the existence of the multipher do not depends on the 
potential V. 

There are many examples of the Chaplygin reducing multiplier for k — 2. Since many 
conditions on the metric and constraints are imposed, until recently there were no nontrivial 
examples of multidimensional systems, appart of several examples for fc = 3, 4 with the prop- 
erty that factor M{q) depends only on one coordinate, that are reducible to a Hamiltonian 
form by the Chaplygin procedure ( gSI EH EHI El ) • 

As an alternative, in the reduction of Chaplygin systems one can use the symplectic (or 
Poisson) framework (see |SDj^ JJ_,_14 ). Such systems can be represented in a Hamilton-like 
form with respect to an nondegenerate (almost-symplectic) 2-form, which however may be 
not closed Namely, let 5 be the Legandre transformation of the semi-basic form S. Then 
one can write H1.9|l as 

^nh{XH,-) ^ dH{-), where fi„h = + S. 

In this framework, the Chaplygin multiplier is a function M such that the form 57 = Mflnh is 
closed. Then, after rescaling Y = X/Af, we obtain the Hamiltonian system Cl{Y,-) ~ dH{-) 
(see |Sni El El E] ) ■ Contrary to the procedure described in Theorem ll.2l here the vector 
field Y has no direct mechanical description. 

Recently, necessary and sufficient conditions for the existence of an invariant measure of 
the reduced system in case when the Lagrangian of the system is of a pure kinetic energy 
type are given in ^1 El • 

1.3 Contents of the Paper 

In section 2 we consider the systems with left-invariant metrics and left -invariant con- 
straint distributions, so called LL systems. The equations of the motion reduce to the 
Euler-Poincare-Suslov equations on the corresponding Lie algebra. Although such equa- 
tions generally are not Hamiltonian, their nice algebraic structure allows us to construct 
various integrable examples with an invariant measure. 

In section 3 we consider a class of LR systems (left-invariant metrics and right-invariant 
constraint distributions), which can be regarded as Chaplygin systems on the principle 
bundle G Q = G/H, H being a Lie subgroup. We show that, in contrast to generic 
Chaplygin systems, the reductions of our LR systems onto the homogeneous space Q always 
possess an invariant measure. Then we study the case G = SO{n), when LR systems 
are multidimensional generalizations of the Veselova problem of a nonholonomic rigid body 
motion, which admit a reduction to the system with an invariant measure on the (co)tangent 
bundle on the unit sphere S"^^. For a special choice of the left-invariant metric on SO{n), 
we prove that under a time reparameterization, the reduced system becomes an integrable 
Hamiltonian system describing a geodesic flow on the unit sphere S"~^. This provides a 
first multidimensional example of a nonholonomic system for which the celebrated Chaplygin 
reducibility theorem is applicable. Lastly, we present an explicit reconstruction of the motion 
on the group SO{n). 
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Finally, in section 4 we present another class of systems on an unimodular Lie group G, 
which always possess a non-trivial invariant measure and which are obtained as modifications 
of a geodesic flow on G with respect to a sum of a left- and a right-invariant metrics, so called 
L-l-R systems. It appears that a nonholonomic LR system on a group G can be obtained 
as a limit case of an appropriate L-l-R system on this group. As an example, we consider a 
nonholonomic mechanical system called the spherical support. 

2 LL Systems 

2.1 Euler— Poincare-Suslov Equations 

In this section we consider nonholonomic systems (G, D) with a left-invariant distributions 
D and a left-invariant Lagrangians I that describes left-invariant metrics on a compact 
connected Lie group G. Let q = TjdG be the Lie algebra of G. In what follows we shall 
identify g and g* by Ada invariant scalar product (•, •), and TG and T*G by bi-invariant 
metric on G. For clearness, we shall use the symbol oj for the elements in g and the symbol 
X for the elements in g* = g. 
Let 

5 = {tj e G, (w, a*) = 0, i = 1, . . . , p} C g 

be the restriction of the left-invariant distribution D to the algebra, for some constant and 
linearly independent vectors a* in g. From the left invariance condition we have Dg = (7 • 0. 
The distribution is nonintegrable if and only if is not a subalgebra. Also, it is sufficient 
to give a Lagrangian at one point of the group, for instance the identity l{g,g) = ^{Iuj,u!), 
OJ = ■ g. Here X : g — > g is a symmetric positive definite (with respect to (•, •)) operator. 
The Hamiltonian in the left-trivialization is given by H{x) = ^{A{x),x), A = I^^. The 
corresponding left-invariant metric will be denoted by dsj. 

Let m be the restriction of the constraint submanifold M to g, that is m = 1(5). Equa- 
tions are G-invariant and reduce to m, 

p p 
X = [x, VH{x)] + A,a* = [x, A{x)] + ^ \^a\ (2.1) 

1=1 i=l 

where Ai are Lagrange multipliers chosen such that x belongs to m = 2^(5), i.e., such that lo = 
A{x) belongs to 0: (-4(a;), a') = 0, i = 1, . . . , p. In other words, the following commutative 
diagram holds 



M - 




^ M 


A 








A 






p* 






m 




-> m. 



where T'* and P* are phase flows of the nonholonomic geodesic flow and the system H2.1f) 
respectively, and A maps g ■ x E TgG to x £ q. 

Following [23, we shall call l|2.1|) the Euler-Poincare-Suslov (EPS) equations, as a gen- 
eralization of the Suslov problem of the nonholonomic rigid body motion (sec the example 
below). 

These equations have a quite different nature in comparison with the Euler-Poincare 
equations x = [x, A{x)]. In particular, as indicated in 41 , in the case of only one constraint 
{a,A{x)) — 0, they have a smooth invariant measure if and only if [a,yl(a)] = /ia. 
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Reconstruction of the motion on the group. In the Hamiltonian case, the integrabil- 
ity of the reduced system imphes generally a non-commutative integrability of the original 
system, namely the phase space is foliated by invariant isotropic tori with quasi-periodic 
dynamic (see |(j(J)'). However there is no such analog in the nonholonomic setting. To re- 
construct the motion (g{t), g{t)) on the whole phase space, we have to solve the kinematic 
equation 

g-\t) ■ g{t) ^ ^(t) ^ A^it)). 

where x{t) are solutions of H2.1|l . i.e., to find all trajectories in A4 that projects to the 
given trajectory x(t) in m. In particular, if x{t) is a relative equilibrium {x{t) = x(to) for 
all t) or if x{t) is a relative periodic orbit {x{t + T) = x{t) for all t), then the invariant 
set A^^{{x{t),t e R}) C M is foliated by invariant tori of maximal dimension rankG or 
rankG + 1, respectively (e.g., see fI7\). 



Multidimensional Suslov problem. The most natural example of LL systems is the 
nonholonomic Suslov problem, which describes the motion of an n-dimensional rigid body 
with a fixed point, that is, the motion on the Lie group SO{n), with certain left-invariant 
nonholonomic constraints. 

For a path g{t) e SO{n), the angular velocity of the body is defined as the left- 
trivialization uj{t) — 9^^ ' 9{t) G so{n). The matrix g G SO{n) maps a coordinate 
system fixed in the body to a coordinate system fixed in the space. Therefore, if ei — 
(eii, . . . , ein)^ ■ • ■ , Cn = (e„i, . . . , e„„)"^ is the orthogonal frame of unit vectors fixed in the 
space and regarded in the moving frame, we have 



Ei^ g- ei, 



En = g ■ e-r. 



where i^i = (1, 0, . . . , 0)^, . . . , i?„ = (0, . . . , 0, 1)^. From the conditions Q = Ei = g ■ Ci 
g ■ Ci, we find that the vectors ei, . . . , e„ satisfy the Poisson equations 



+ 



-UJCi 



1, 



(2.3) 



The left- invariant metric on SO{n) is given by non-degenerate inertia operator X : 
so{n) — > so{n). Then the Lagrangian of the free motion of the body reads I — ^{Iuj,uj), 
where now (•, •) denotes the Killing metric on so{n), {X,Y) = — -ifr (XY), X,Y £ so{n). 
For a "physical" rigid body, luj has the form luj -f loI , where / is a symmetric n x n matrix 
called mass tensor (see |24) ) . However, since we are interested mainly in nonholonomic 
geodesic flows, we shall consider other inertia operators as well. 

Recall that in the three-dimensional case, the Suslov problem describes the motion of a 
rigid body with the constraint: the projection of the angular velocity to a vector fixed in the 
moving frame (for example £^3) is equal to zero |52[ 1^ . In other words, only infinitesimal 
rotations in the planes sj>a'n{Ei, E2) and span(i?i, £"3) are allowed. Hence, it is natural to 
define its n-dimensional analog as follows: only infinitesimal rotations in the fixed 2-planes 
spanned by (i?i, £^2), • ■ • , {Ei, En) (i.e., in the planes containing the vector Ei) are allowed. 
Following gH, one can relax these constraints by assuming that the angular velocity matrix 
has the following structure 



/ 



UJir 



OJln\ 



-LOlr 



• • • LOrn 

o 
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where O is zero (n — r) x (n — r) matrix. This imphes the left -invariant constraints 

(w, EiAEj) ^0, r+l<i<j <n. (2.4) 
As a result, the Suslov problem is described by the EPS equations 

^ (lio) = [luj, uj]+ J2 Ep^^i^ (2-5) 

r<p<q<n 

together with Poisson equations (|2.3|) . Here the components of the vectors ei, . . . ,e„ play 
the role of redundant coordinates on SO{n). 

Various integrable cases of the Suslov problem with additional potential fields and their 
multidimensional generalization are given in [^OHIISEI and respectively. 

2.2 Some Integrable Cases of EPS Equations 

EPS equations on symmetric pairs. Let f) be the subspace of the algebra g spanned 
hy a\ i=l,...,p. 

Consider the case when the tensor A preserves the orthogonal decomposition g = f) + 0, 
i.e., A = Atj + Ao, where Af, : f) ^ (), Ac : t) ^ D are positive definite operators. Then 
m = /(O) = 5, and we can write (|2.1|l in the following way 

X ^ [x,Azi{x)]^, xGd, (2.6) 

where denotes the orthogonal projection of ^ G g to the subspace (with respect to (•, •)). 

Equation H2.6|) preserve the standard measure on 5. Also the Hamiltonian function 
H{x) = ■^{x, Ad{x)) and the invariant F{x) = {x, x) are always first integrals of the system. 
Therefore, by the Jacobi theorem the equation (|2.6|l is always integrable for dimO < 4. 

Remark 2.1. Note that, in general, the invariant F{x) — {x,x) is not the integral of (|2.1|l . 
although it is always an integral of non-constrained system. Namely, a first integral f{x) of 
the Euler-Poincare equations x = [x, A{x)] is the integral of (|2.1|) if and only if the following 
condition holds 

^A,(V/(x),a*)Uen = 0. (2.7) 

i 

In our case WF{x) = 2x, a; e m = is orthogonal to t) and therefore the invariant F{x) 
remains to be an integral. 

Example 2.1. Let 6 be a subalgebra of g and tr the orthogonal complement of i. Suppose 
that (g, £) is a symmetric pair, i.e., the following conditions are satisfied: 

[ft),tr]Ct, [6,t]Ct, [t,tt3]CtT3. 

Then, in the special case c) = fo, we have [O, t)]o ~ 0. Therefore all the solutions of equations 
1)2 .611 are constants. As a result, the solution of the original system on G (nonholonomic 
geodesic lines of the metric dsj) is given by the motion along one-parameter subgroups, 

g{t) = gacxp{t^), ^ G 5. 

This simplest situation occurs for the multidimensional Suslov equations (|2.5I) with r = n—1 
and Juj = lu + tol, where / = diag (/i, . . . , /„). Then f) = so{n — 1), (so(n), so{n — I)) is 
a symmetric pair, and X preserves the decomposition so{n) = + so{n — 1). Hence the 
solutions ijj(t) are just constants. 
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Motivated by the above observation and by another integrable case of the multidimen- 
sional Suslov problem (see below), let us assume that there is a chain of subalgebras 

where (g, t) is a symmetric pair, and consider the adjoint representation of 1 onto the linear 
space xv: rj G I ^-^ [rj, ■] G End{'m). With respect to this representation, decompose to into 
irreducible subspaces ro — roo + Wi + • • • + tOm, tuo being the subspace with the trivial 
representation. Next, assume that has the form 

c) = u+c)o+fiH \-'0g, u^dnl, dk='OntVk. (2.8) 

with dimu > 1, g < m. Suppose also that = s ■ Id^, s G M and that the operator Au 
preserves the decomposition (I2.8|l . that is, Aj, = + Ao + ■ ■ ■ + Ag. Let Xk denote the 
orthogonal projection of a; to 0^, k — 0, . . . , g. Then the equation (|2.6() reads 

■^ixu+ xo-\ \- Xg) ^ [xu + xo ^ \- Xg, sxu + Ao{xo) H h Ag{Xg)]o. 

In view of conditions [u, Ofcjg C dk, [fj, f j]o C u, the above system splits into g + 2 equations 
io = 0, 

iu = N H \- Xg,Ao{xo) H h Ag{xg)]^, 

Xk ^ [xu,Bk{xk)]i,^, /c = l,...,5, 

where Bk = Ak — s ■ Idk ■ 

Thus, apart from H{x) = ^{Ai3{x),x) and F{x) — {x,x), the system (|2.6|l has a set of 
the first integrals given by the projection of x to do, Fq{x) — xq and the functions 

Fk{x) ^ {Bk{xk),Xk), k^l,...,g. 

In this case the following theorem holds (see 

Theorem 2.1. 1). If the operators Bk are positive definite, then invariant varieties 

Mc ^ {x G D \ xo ^ co,Fi(x) = ci, . . .,Fn{x) = Cg, F{x) = Cg+i}, 

Ci, . . . , Cg+i being constants of motion, are diffeomorphic to the product of spheres 
^.dimoi-i X ... X 5"iimDg-i ^ S^^™""!, provided that Cg+i satisfies inequality 

a 

cg+i > |coP+^^, bk^ mm {Bk{^k),^k). (2.9) 
fr[bk I?. 1=1 

2) If dim dk < 2, k — \, . . . ,g, dimu — 1, and all the constants Ci are nonzero, then Mc 
is diffeomorphic to the disjoint union of two g -dimensional tori with a quasi-periodic 
dynamic of the form ipi = Qi/^{ip), i = 1, . . . ,g. 

The Fedorov— Kozlov case. The above construction applied to the symmetric pair (g, t) — 
{so{n),so{2) X so{n — 2)) gives the Fedorov-Kozlov integrable case of the multidimensional 
Suslov problem (5^. As above, let in be the orthogonal complement of t: 

so{n)=(^';^^^} ,0(^-2))- 
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We take u = I = so(2) = span{£'i A £'2}, 3 = u + ro, i.e., the constraint are given by relations 
(E3I) with r = 2. Then 



Oi = iTi = span{i;i A S3, A E3}, dn-2 = n)n-2 = span{£'i A En,E2 A En}. 

In the Suslov problem there is a natural choice of the inverse inertia operator A = T^^ 
which preserves the decomposition 5 — so{2) + Oi + • • • + fn-2- Namely, we take the left- 
invariant metric on so{n) determined by the kinetic energy of the multidimensional rigid 
body, 

A:E.^A Ej ^ Y~~T^' ^ ^3 ■ 
li + Ij 

li Ii > I2 > I3 > ■ ■ ■ > In, then, under condition ()2.9|) . the integrals Fi,...,F„_2 are 
positive definite and the invariant submanifolds 

{Fi = Ci, . . . , F„-2 = Cn-2, F = c„_i} 

are union on two disjoint (n — 2)-dimensional tori. Moreover, as shown in )24| . the motion 
on the tori is straight-line but not uniform and in appropriate angle coordinates (fi it is 
described by equations 



O Kh - Ii+2)(l2 - Ii+2) . , „ ,'0 11^ 

t^uiv) V (-^1 + h+2){l2 + I1+2) 

For this integrable case the reconstruction problem was studied in |2H1 • As follows from 
(|2.11|l . if the trajectories are periodic on one torus, they are periodic on the rest of the tori. 
Then the trajectories {g{t),g(t)) which correspond to the given periodic trajectory uj{t) are 
quasi-periodic (e.g., see pTjl. 

According to (SS], in the opposite case, if for some constants c > and 7 > n — 3, the 
frequencies satisfy Diophantine conditions 

\l + V^{k,n)\>c/\k\'', ; = 0,1,2, foraU fc G Z"-^ 

and the value of the integrals F (or F and Fi) are dominant with respect to those of other 
integrals (ci/c„_i ~ e, i 7^ 71 — 1 or Ci/ci, Ci/cn-i ~ e, i ^ 1, n — 1), then the dynamics on 
the whole phase space can be approximated by quasi-periodic dynamics on the time interval 
of length ~ exp(l/e). 



The Suslov problem on so(4). Now we concentrate on the integrable case when is 
not an eigcnspace of A. Let g = so(4). Then { = span{i?i A E2,E3 A E4} is a Cartan 
subalgebra. As above, take the inertia operator in the form (|2.1U|I . which implies A{i) = t. 
Further, take a — aiEi A E2 + a2i?3 A E4, aia2 ^ and the constraint (a, A{x)) = 0. Then 
[a, -4(a)] — and the Euler-Poincare-Suslov equations 

X — [x, A{x)] + Xa 

preserve the standard measure on m = {a; | {a,A{x))} = (see j^). Next, one can 
always chose a linear combination of quadratic invariants on so(4), ci/i -I- C2/2, such that 
the condition H2.7fl holds. Thus our system on five-dimensional space m has the integrals 
H = ^{A{x), x) and Fi = ci/i -f C2/2. For the integrability one needs one more independent 
integral. It can be taken in the form of a quadratic function on the orthogonal complement 
oft 

This approach is a special case of the method of construction of integrable EPS equations 
on six-dimensional unimodular Lie algebras given in |3l)|. 
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Chains of subalgebras. Suppose there is a chain of subalgebras 

00 C fli C • • • C s„ = 0. 

Let Qi = 0j_i + roi be the corresponding orthogonal decompositions. Then fli = 0o + toi + 
■ ■ ■ + tOi. Following [7], consider A of the form: 

A = Ao + si ■ Idm -\ \-Sn-Idn,^, Si > 0, i = l,...,n, (2-12) 

where is a symmetric positive operator defined in the subalgebra go- Suppose that has 
orthogonal decomposition 

O^Oo+fiH hOn, Ofc = {tjfc e tTfe, (afc,Wfc) = 0, i = 1, . . . (2.13) 

Then D/;, A; > are invariant subspaces of A. By Xk denote the orthogonal projection of x 
to Ofc, fc > 0; and by xq denote the orthogonal projection to go- 
Now we can formulate the following theorem (see ISlp. 

Theorem 2.2. The Euler-Poincare-Suslov equations 1^2.1]) . with D and operator A{x) of 
the form \2.1!^) and \2.1^) . are equivalent to the Euler-Poincare-Suslov equations on the 
Lie subalgebra go-' 

Po 

io = [a;o,-4o(a;o)] +^Mi«o> (2-14) 

1=1 

(A(a;o),ao) =0, i = l,-..,po, 
together with a chain of linear differential equations on the subspaces 5^; 

Xk = [xk.Aoixo) - SkXo + (si - Sk)xi H h (sfc-i - Sk)xk-i]k- (2-15) 

If the Euler-Poincare-Suslov equations (|2.14l) on go are solvable, then the integration of 
original equations (|2.1() reduces to consequitive integration of the chain of linear dynamical 
systems (|2.15ll for fc > 0. In the most simplest case the solutions of (|2.14ll are constants. 
Then the components of the vector xi satisfy a system of linear equations with constant 
coefficients, hence they are elementary functions of the time t. This happens if ^o -^c^go 
or if go is a commutative subalgebra. In particular, if dim Z?o = 0, then we have a;o = 0. In 
this case ii = and X2 is given by elementary functions of t. 

2.3 Hamiltonian Flows 

In some cases, the nonholonomic geodesic flow (|1.3ll on Ai can be obtained as a restriction of 
a Hamiltonian flow on T*Q to the invariant submanifold M. In most examples this happen 
when the Lagrange multipliers in (|1.3|l vanish, i.e., when M is an invariant submanifold of 
the unconstrained geodesic flow. 

There are also cases of nonzero Lagrange multipliers. This means that M is the invariant 
submanifold of some other Hamiltonian system. In particular, in Example l2.1l one can take 
a geodesic flow of a bi-invariant metric. Note that the Lagrange multipliers, in general, are 
different from zero (X]f=i •^i'^* ~ "[2^7-^0(2^)]) ^-i^d is not an invariant submanifold of the 
unconstrained geodesic flow of the left -invariant metric dsj. 

Below we concentrate on the first case (zero Lagrange multipliers). Suppose that the 
orthogonal complement f) of J) is the Lie algebra 6 of a Lie subgroup and that the operator 
A also preserves orthogonal decomposition g — f) + f), i.e., A — Atj +Atj- Then m = 1(5) = c) 
and EPS equations take the form 1)2. 6|l . 
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Further, suppose that {x,A!3{x)) is an invariant of the adjoint action of ^ on 0, 

{[t),x],Ai{x)) — 0, or, equivalently, [x, ylo(a;)]f| =0, a: e 0. 

Then one can easily check that 5 is the invariant subspace of the Euler equations 

X = [x, VH{x)] = [x, A{x)], xeQ. (2.16) 

Therefore the Lagrange multiphers vanish. Moreover, one can consider the new Hamiltonian 
function H*{x) = ^{x, A-!}{x)) on g and the Euler equations 

x=[x,'^H*{x)]^[x,At,{x)], xeQ. (2.17) 

In both cases, the restriction of the systems to coincides with the Euler-Poincare-Suslov 
equations (|2.6|) . For example, after projection to f) and the system (I2.17f) becomes 

ij, = [xx, , At,{xj))] + [xt,At,{xt))], if, = 0. 

Let ft., h* : T*G ^ R be the functions obtained by left translations from H and H*. 
While h is the Hamiltonian of the geodesic flow of the left-invariant metric dsj, the function 
h* is degenerate in momenta and has another geometric meaning. 

Suppose that d generates the Lie algebra g by commutations. Then, by the Chow- 
Rashevski theorem, any two points on G can be joined by a piecewise smooth admissible 
curve g{t). Locally shortest admissible curves are called sub-Riemannian geodesic lines of 
the sub-Riemannian metric obtained by restriction of the given left-invariant metric dsj to 
D. The Hamiltonian flow of h* on T*G is a sub-Riemannian geodesic flow. In other words, 
the projection of the flow to G give us sub-Riemannian geodesic lines (for more details see 
|51[l53j ). Such systems are also known as vaconomic systems |2- 

We summarize previous considerations in the following proposition (see |31p. 

Proposition 2.3. Suppose that H\^ is an invariant of 9j-adjoint action and thatt) generates 
the algebra g by commutations. Then on the constrained submanifold A4 the following three 
different problems have the same flow: the nonholonomic geodesic flow, the geodesic flow 
with Hamiltonian h and the sub-Riemannian geodesic flow with Hamiltonian h* . 

An example on the Lie group SU(n). Let us illustrate how the special case of the 
construction given in the Theorem 12.21 produces a nonholonomic geodesic flow with the 
above property. Namely, consider the chain of subalgebras 

su{2) C su(3) C • • • C su{n) 

given by the natural matrix embedding. Let su{2 + i) — su{2) -f roi be the orthogonal 
decompositions and let A has the form 

A^ S^Idsu{2) + Sl ■ Idn^ ^ V Sn-2 ■ Idxr,,,_2, Si > 0, i = 0, . . . , 71 - 2 . (2.18) 

We take i) to be the orthogonal complement to the Lie algebra of the maximal torus T"^^ C 
SU{n) consisting of diagonal matrices. Then the Hamiltonian H = ^{x,Ax) will be an 
invariant of the adjoint action of T"^^ on su{n) (see, e.g., ^Ul) and D will generate so{n) 
by commutations. Thus the system satisfies the conditions of ProDOsition l2.3l 

Furthermore the system is integrable and can be considered as a Chaplygin system as 
well. Namely, let h be the corresponding left invariant Hamiltonian function on T*SU{n). 
Since H is adjoint T"^^ -invariant, we have that h is also right T"^^ -invariant function. 
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Thus, the group T" ^ acts on Riemannian manifold {SU{n),dsj) by isometries. By sub- 
mersion, the metric dsj, induces the S'[/(n)-invariant metric dsj^^^^ on the flag manifold 
Fn = SU{n)/r-^ 

T"-i — > SUin) 

i TT . (2.19) 

F„ = 5(7(n)/T"-i 

Note that the horizontal spaces of the submersion coinside with those of the distribution D. 
In other words, we can also consider {SU{n),ds'^,D,T^^^) as an example of a Chaplygin 
system such that the right hand side of H1.7|l is equal to zero (see Remark ll. If) . 

The geodesic flow of the metric ds^^^f^ on F„ is completely integrable (see [TOj)- To 
describe the motion on the whole phase space D one must solve the reconstruction problem. 
Since the group T"^^ is Abelian, this can easily be done by quadratures (e.g., see |42)'). 

3 LR Systems 

3.1 LR Systems as Generalized Chaplygin Systems 

Following j54[ I55| , one defines an LR system on a compact Lie group G as a nonholonomic 
Lagrangian system (G, I, D) where I is a left-invariant Lagrangian and Z? is a right-invariant 
distribution on TG. As in LL systems, the Lagrangian is defined by a left-invariant metric 
dsj, l{g,g) = \ {Iu,uj), uj ^ g''^ ■ g. 

The distribution D is determined by its restriction d to the Lie algebra, 

Dg=d-g^g- {g-' ■d-g)ci TgG. 

Let f) = spanjo"'^, . . . , a^} be the orthogonal complement of d with respect to (•, •). Then the 
right-invariant constraints can be written as 

ujeg^^-(}-g, or (w, 5"^ • a' • g) = 0, i = l,...,p, 
or, equivalently, {a\A{x))—0, = g^^ ■ ■ g. (3.1) 

Equations H1.3|) in the left trivialization take the form 

p 

Aix)]+J2>^^c,\ (3.2) 

i=l 

g = g-uj = g- A{x). (3.3) 

Here the Lagrange multipliers Xi are determinated by differentiating the constraints. 

The system (13.2(1 . ((3.3(1 is actually defined on the whole phase space TG and has first 
integrals 

hi9,x) = {A{x),g^'^ - a' ■ g), i^l,...,p. 

Then the nonholonomic geodesic flow is just the restriction of 1(3. 2() . ((3.3(1 onto the invariant 
submanifold M = {{g,x) \ fi — 0, i — 1, . . . , p}. 

Instead of 1(3. 2() . ((3.3(1 . one can consider the following closed system on the direct product 
0^+^ in the variables {x, , . . . , a^}, 

p 

x=[x,A{x)]+J2^^'^'' (3-4) 

1=1 

a' = [a\Aix)], i^l,...,p, (3.5) 
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where the multipliers are determined from the conditions ■^{a'^,A{x)) = 0. 

Equations (|3.5|) imply that a^{t) belongs to the adjoint orbit OG(a*(to))- Then, if 
{x{t),a^{t), . . . ,aP{t)) is a solution of (|3.4|) . (|3.5|) and g(t) is a solution of the kinemati- 
cal equation (|3.3|) (with appropriate initial conditions), we conclude that {g{t), x{t)) is a 
solution of the system H3.2|l . (|3.3|l . 

One of remarkable properties of LR systems is the existence of an invariant measure, 
which puts them rather close to Hamiltonian systems. Veselov and Veselova j55| proved 
that the system (|3.4ll . (|3.5(l has an invariant measure with density 

y/det{{Aia'),ao)). (3.6) 

This implies that the original system H3.2|l . (|3.3(l on TG also has an invariant measure of 
the form fJ.{g) ■ da, where da is the canonical volume form on TG and 

^l{g) = det {{A{g-^ ■ ■ g), g-^ ■ -g)). 

In particular, our nonholonomic geodesic flow on A4 also has an invariant measure described 
in the following way. Let d S be a volume form on M . Then 

da dfiA--- AdfpAdY., {g,x)eM (3.7) 

for some positive function 6. Next, let £ be the Lie derivative with respect to the flow (|3.2|l . 
(13.3(1 . Since the functions fs are first integrals, we have £dfs = 0, s = 1, . . . , p. As a result, 
from the condition C{fj, da) = and (|3.7|l we obtain d/i A • • • A dfpC{fi9 dE) = 0. Hence, the 
restriction of the flow onto A4 has the invariant measure fj,9 dYi. 

Reduction. Now, let the linear subspace f) be the Lie algebra of a subgroup H C G. 
Then the Lagrangian l{g,g) and the right-invariant distribution D are also invariant with 
respect to the left iJ-action. Consider homogeneous space Q = H\G of cossets {Hg}. The 
distribution D can be seen as a principal connection of the principal bundle 

H — > G 

Q = H\G 

As a result, the LR system (G, Z, H) can naturally be regarded as a generalized Chap- 
lygin system. In order to write the reduced system on Q in a simple form, we identify g 
and g* by the Adc-invariant scalar product (•,•), and the spaces TQ, T*Q by the normal 
metric, induced by the bi-invariant metric on G. Next, consider the moment mappings: 

(b:TG^T*G^2, <i> ■ TQ - T*Q ^ 2, 

of the natural right actions of G on T*G and T*Q, respectively. We have (l){g) = uj = g^^ ■ g 
and the map $ can be considered as a restriction of to D. 

The reduced Lagrangian is, by definition, the constrained Lagrangian 

considered on the orbit space H\D = T{H\G). It follows that the reduced Lagrangian is 
simply given by 

Liq,q) = i(X$(g,g),$(g,(z)), 

where q — Tr{g) are local coordinates on Q (which may be redundant). This is a Lagrangian 
of the geodesic flow of metric which we shall denote by dsj ^. 

By using equations (|1.7|) one can prove the following proposition (see |24| ) . which is a 
special case of the general nonholonomic reduction procedure described in [361 |Sj . 
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Proposition 3.1. The reduced Lagrange-d'Alambert equation describing the motion of the 
LR system {G,l,D) has the form 

^-j^^,^^^{I<fiqA),W,-r^,mq,q),Hq,0]), (3-8) 

for all virtual displacements ^ e TqQ, where pr^-^g : g — )■ g~^i)g is the orthogonal projec- 
tion, and q = 7r(g). 

In addition, it appears that tlie reduced LR system 13.8|l also possesses an invariant mea- 
sure (note that a generic Chaplygin system does not have this property, see Namely, 
the following general statement holds (e.g., see |26|). 

Lemma 3.2. Suppose there is a compact group acting freely on a manifold N with local 
coordinates z and there is a ^-invariant dynamical system z — Z{z) on N. If this system 
has an invariant measure (which is not necessary <3 -invariant), then the reduced system on 
the quotient manifold N / also has an invariant measure. 

3.2 Veselova Problem, an Integrable Geodesic Flow on the Sphere 
and the Neumann Problem 

Veselova problem. The most descriptive illustration of an LR system is the Veselova 
problem on the motion of a rigid body about a fixed point under the action of nonholonomic 
constraint 

(f^,7) = 0, (3.9) 

where e M"^ is the angular velocity vector, 7 e M"^ is a unit vector, which is fixed in a 
space frame, and (, ) denotes the scalar product in 'E? [SJ. Geometrically this condition 
means that the projection of the angular velocity of the body to a fixed vector must equal 
zero. 

The equations of motion in the moving frame have the form 

Xf2 = Xf7 X SI + A7, 7==7xf7, (3.10) 

where I is the inertia tensor of the rigid body, x denotes the vector product in R'^, and A 
is a Lagrange multiplier chosen such that Vl{t) satisfies the above constraint, 

(I ^7,7) 

The Veselova system (I3.9II . (|3.10|l is an LR system on the Lie group S'0(3), which is the 
configuration space of the rigid body motion. After identification of Lie algebras (R^ , x ) 
and (so(3), [•,•]), the operator I induces the left-invariant metric ds']. The angular velocity 
correspond to 17 = g~^g, the velocity in the left trivialization TS0{3) ^ 5*0(3) x so(3). The 
vector fixed in the space corresponds to the right-invariant vector field Jg = g ■ ig~^ ' ' 9) G 
TgSO{3), a € so(3), and the nonholonomic constraint (|3.9(l has the form {g~^ ■ a - g,fl) =0. 
Once can check that the closed system (|3.1UII . (|3.11ll has invariant measure with density 
(7,X~-'^7), as predicted by formula (|3.()|l . Note that integrable potential perturbations of the 
Veselova system can be found in ,54| |26j . 

Multidimensional Veselova problem. Now we proceed to a n-dimensional generaliza- 
tion of the Veselova system, describing the motion on the Lie group SO{n) with certain 
right-invariant nonholonomic constraints. 
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Let, as above, ei, . . . , e„ be unit vectors that form a fixed orthogonal frame in the space 
M". Then, similarly to the generalized Suslov problem in Section 2, we define n-dimensional 
analog of (|3.9|l as follows: only infinitesimal rotations in the fixed 2-planes spanned by 
(ei, 62), . . . , (ei, e„) are allowed. This implies the constraints 

{uj, A Cj) = 0, 2<i <j < n. (3.12) 

Equivalently, consider the right-invariant distribution D on TSO{n) whose restriction to 
the algebra so{n) is given by 5 = span{£'j h Ek, k — I, . . . ,r, j = I, . . . , n}, where Ei A Ej 
form the basis in so{n). Since Ci A ej = g^^ ■ Ei A Ej ■ g, we have that constraints are 
UJ = g~^ - Vi ■ g — span{ei A ei, 2 < i < n}. 

Remark 3.1. As for the multidimensional Suslov problem, the constraints (|3.12ll can be 
relaxed. However, in this case, the existence of the integrable LR system is still not known. 
That is why we keep using the above constraints (see Theorem 13. 4|l . 

The LR system can be described by the Euler-Poincare equations 1)3.21 13.31 13.4|l on the 
space so{n) x SO{n) with indefinite multipliers Apq, 

2<p<q<n 

ei+ujei=0, i = l,...,n. (3.13) 

Here, as above, the components of ei,...,e„ play the role of redundant coordinates on 
SOin). 

Reduction. The orthogonal complement f) of is a Lie algebra, namely 

t) = span{i?p A Eq, 2 < p < q < n} ^ so{n — 1). 

Therefore, the Veselova system can be treated as a generalized Chaplygin system on the 
principal bundle 

50(71-1) — > SO{n) 

i , (3.14) 

S*"-! = SO{n - l)\SO{n) 

where 5"^^ is the n-dimensional sphere, realized as the unit sphere in M", 

5"-i = {ge M"-\ ql + ---+ql = l}, 

where we set q = d. The moment map is then uj = ^{q, q) — qAq. Thus, for solution ei{t), 
U!{t) = ei{t) A ei{t) of l|3.13(l . q{t) = ei{t) is a motion of a reduced system on the sphere 

The invariant measure. It appears that for some special inertia tensors, many of the 
calculations takes an especially simple form. Suppose that the operator X is defined by a 
diagonal matrix A = diag(j4i, . . . , An) in the following way 



I(E,AE,)^^^E,AE,. (3.15) 



Notice that for n = 3 this corresponds to the well known three-dimensional vector formula 
I{x X y) = {detA)~^ Ax x Ay, A = I^^. 
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Under the condition H3.15|l the reduced Lagrangian L{q,q) and the right hand side of 
the Lagrange-d'Alambert equation (|3.8|) take the form 

^ [iAq,q)iAq,q)^{Aq,qf], (3.16) 



2det^ 

q), prg-i^g[$(g, q),<^{q, 0]) = ^^(^9 ^ Ag, pr^-i^^ ^ A 

Here we used relation pr^-i^g ^ A g = ^ A g for any admissible vector ^ = (^i, . . . , ^„)-^ G 
TqS"'~^. Below we shall keep using the redundant coordinates qi and velocities qi, in which 
the Lagrange equations have the form 

dL d dL /o 
a 17"a^ = + i = l,...,n, (3.18) 

dgi di dqi 

0* 1 1 

where A is a Lagrange multiplier. 

Now we want to represent the reduced LR system on T*S'^~^ as a restriction of a system 
on the Euclidean space K.^" = {q,p}. Note that L{q,q) is degenerate in the redundant 
velocities q, hence they cannot be expressed uniquely in terms of the redundant moments 

P. = g ^ ^il^M)A.q. - ^i^.Aq)Aq^. (3.19) 

In this case one can apply the Dirac formalism for Hamiltonian systems with constraints in 
the phase space (see, e.g., ^1I21ES1)- Namely, from ()3.19(l we find that {q,p) — 0, hence the 
cotangent bundle r*5'"~i is realized as a subvariety of M^" — {q,p) defined by constraints 

01 = {q, q) = 1, 02 = {q,p) = 0. 
Under these conditions, relations (|3.19l) can be uniquely inverted to yield 

det^l 



(g, Aq) 



[A-'p-{p,A-\)q] . (3.20) 



On the other hand, we note that dL/dqi = tt^. Then, from H3.18|l we obtain p — —Aq and, 
from the condition {q,p) + {q,P) — 0, 

p=-Kq, A = detA^^i^^:^^^)^i#(^l^. (3.21) 

{q,Aq) 

The system (I3.20f) . ()3.21|l on T*S"^^ coincides with the restriction of the following 
system on R^" — {q,p} 

q^^{qi,H}^, Pi = {pi,H}:,-TTi, 

1 (p A~^p) 
TT,{q,p) ^ Tr,{q,q{q,p)), H ^ - detA , 

which is quasi-Hamiltonian with respect to the following Dirac bracket on R^" 

{01, 02} 
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{•, •} being the standard Poisson bracket on R^". This system has exphcit vector form 

detA 



(q, M) 



, , , ip,A-^p){q,q)- {p, q){q,A-^p) 
raw ' 



(3.22) 



The bracket {•, is degenerate and possesses Casimir functions 0i, (^2 specified above. 

Now, we can find the exphcit form of the invariant measure of the reduced system. From 
(IT^ we find 

^^{dqr dpj ^" > {q,q)iq,Aq) ' 

which, in view of H3.19|l . takes the form {n — 2){q, Aq)/ {q, Aq). Hence the extended system 
()3.22|l possesses an invariant measure 

J = {Aq, g)-("-2)/2 dpi A A • • • A dpn A dg„. 

Next, at points of T*S^^^^, the standard volume form in M^" can be represented as 

dpi A dqi A ■ ■ ■ A dpn A dqn = w""^ A A (i(f>2, 

where w is the restriction of the standard symplectic form dpi A dqi + • • • + dpn A dqn onto 
rp*gn-i g^^^ $i,$2 are certain functions of the Casimir functions (f>i,4>2- Since the latter 
are invariants of the vector field V{p, q) given by (|3.22l) . the Lie derivatives Cyd^i, Cvd^2 
equal zero. Then, since LvJ — 0, we conclude that on T*S'"^^, 

As a result, we arrive at the following theorem. 

Theorem 3.3. The reduced LR system \S.2(\ \!i.21\) on T*S"'~^ possesses an invariant 
measure 

f{q) = {Aq, g)-("^2)/2 ^ ^ ^„-l 

where a is the canonical volume 2(n — \)-form on T*S"^^. 

Chaplygin reducing multiplier. As follows from Theorem 13.31 item 1) of Theorem 
11.21 and the fact that the dimension of the reduced configuration manifold equals 71 — 1, if 
our reduced LR system on T*S"~^ were transformable to a Hamiltonian form by a time 
reparameterization, then the corresponding reducing multiplier J\f should be proportional 
to l/^{q,Aq). 

Although Chaplygin's reducibility theorem does not admit a straightforward multidi- 
mensional generalization, i.e., item 1) of Theorem 11.21 cannot be inverted, remarkably, for 
our reduced LR system on T*S''-^ the inverse statement becomes applicable (see |2(i|'l. 



Theorem 3.4. 1). Under the time substitution dr — det A/ {Aq, q) dt and an appro- 
priate change of momenta, the reduced LR system or l\y.i!(J\) . 1^3. 2 1\) becomes a 
Hamiltonian system describing a geodesic flow on S"~^ with the Lagrangian 



L*{q,dq/dT)^\{q,Aqr^ 



(3.23) 
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2). For Ai < A2 < ■ ■ ■ < A„ the latter system is algebraic completely integrable for any 
dimension n. In the spheroconic coordinates Ai, . . . , A„_i on S""^^ such that 

2 _ (/i -Ai) ■■■(/, -A„-i) _ 1 

the Lagrangian L*{q,dq/dT) takes the Stdckel form 

1 n.^fc(Afe-A,) /d 

^hi (Afe - /i) • • • (Afc - /„) Afc 
and the evolution of Afe is described by the Abel-Jacobi quadratures 

Xl~^dXi \^J[dX„~l r— /o oc;^ 

, , H \ , = 4.n-i V 2/1 dr, (3.25 

2/R(A0 2Vi?(A„_i) 

A: = 1, • • • , n — 1, 

i?(A) = -(A - /i) • • • (A - /„)A(A - C2) • • • (A - c„_i), (3.26) 

h = L* being the energy constant and C2, • • • ,c„_i fcem;? other constants of motion 
(we set ci = 0). For generic values of these constants the corresponding invariant 
manifolds are (n — I) -dimensional tori. 

The item 1) of Theorem 13.41 is based on the relation between the reduced LR system to 
the celebrated Neumann system (see Theorem 13.51 below*) . 

Namely, consider the iso-energy submanifold £h = {L{q, q) — h} C TS"^^^ of the reduced 
Veselova system ()3.18|) and introduce another new time ti by formula 



dTi = J det A— ^— dt. (3.27) 

V iM, q) 

Theorem 3.5. Under the time substitution \3.21^ , the solutions qit) of the reduced mul- 
tidimensional Veselova system on S"^^ lying on the £h transforms to the solution of the 
integrable Neumann problem with the potential U{q) = ^{A^^q,q), 



d^ 1 

corresponding to zero value of the integral 

2 



-^, = --, + A, (3^28) 



and vise versa. 



For n = 3, Theorem 13. 51 is proved by Veselov and Veselova [HHl- The proof for arbitrary 
dimensions is given in pfij . 



3.3 Reconstructed Motion on D 

Now we consider the integrability of the original (unreduced) LR system on the right- 
invariant distribution D C TSO{n), which is specified by constraints H3.12|l and the left- 
invariant metric given by H3.15|) . The relation between the reduced LR system and the 
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Neumann system described by Theorem l3.5l appears to be useful to reconstruct the motion 
on D exactly. For this purpose we also shall make use of the correspondence between the 
Neumann system and the geodesic flow on a quadric (see Knorrer j35p. Namely, consider a 
family of (n — l)-dimensional confocal quadrics in R", 

g(a) = + ••• + — ^ = -iV "GR. (3.30) 

[a - Ai a - An ) 

Theorem 3.6. (|SS1)- Let X{s) be a geod esic on the quadric Q{0), s being a natural 
parameter. Then under the change of time 



{dX/ds,A-^dX/ds) ^ 
^^^V {X,A-^X) ^^-^^^ 

the unit normal vector q{Ti) = is a solution to the Neumann system 

corresponding to zero value of the integral 2y\} and vise versa. 



It is well known that the problem of geodesies on a quadric Q{Q) is completely integrable, 
and qualitative behavior of the geodesies is described by the remarkable Chasles theorem 
(see e.g., 35, 43^): the tangent line Is = {^{s) + crdX/ds | ct G M} of a geodesic X{s) on 
Q(0) is also tangent to a fixed set of confocal quadrics Q{a2), . ■ . ,(5(q:„-i) C R", where 
a2, ■ ■ • , Oin-i are parameters playing the role of constants of motion (we set ai = 0). Now 
let rifc be the normal vector of the quadric Q{otk) at the touching point = tC]Q{ak). Then 
another classical theorem of geometry says that the normal vectors rii, . . . , n„_i, together 
with the unit tangent vector 7 = dX/ds, form an orthogonal basis in R". 

On the other hand, in |33|, Moser made the following observation. 

Proposition 3.7. 1 ). Let x be the position vector of a point on the line £s, which is tangent 
to geodesic X{s). Then in the new parameterization si such that ds — —(X,A^'^X) dsi 
the evolution of the line is described by the Lax equations in n x n matrix form 



-£=[£,6], C^n-^iA-x(g>x)n^, (3.32) 



d 
dsi 

B = A-^x®A-^j~A-^j®A-^x, (3.33) 



where H-^ = Id — (7, 7) ""^7 7 is the projection onto the orthogonal complement of j 
m R". 

2). The conserved eigenvalues of C are given by the parameters ai — 0, 02, ■ . ■ , ctn-i of the 
confocal quadrics and by an extra zero. The corresponding eigenvectors are parallel to 
the normal vectors xxi — q, . . . , n„_i, and to 7. 

Now we are ready to describe generic solutions of the original LR system on D C TSO{n). 
Let g(ri) be the solution of the Neumann system H3.28|l with Fo(g,q') = 0, which is associ- 
ated to a solution {q{t),p{t)) of the reduced LR system as described by Theorem 13.51 Let 

X = (g, Aq)-'^^Aqis), m = g(s), . . . , n„_i(s), 7(5) = ^ (3.34) 

be the corresponding geodesic on Q{0) in the new parameterization s given by H3.31|l and 
the unit eigenvectors of £. Also, according to H3.27|l and H3.31|l we can treat s as a known 
functions of the original time t. Then we have the following reconstruction theorem (see 

m)- 
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Theorem 3.8. A solution {g{t), g{t)) of the original LR system on the distribution D is 
given by the momentum map uj{t) = q A q and the orthogonal frame formed by the unit 
vectors 

ei = q{t), £2 = n2{t), . . . , e„_i = n„_i(t), e„ j{t). 

The other solutions {g{t), g{t)) that are projected onto the same trajectory {q{t),p{t)) have 
the same uj,ei, while the rest of the frame is obtained by the orthogonal transformations, 

(e^it) ■ ■ ■ e„(0) = (n^it) ■ ■ ■ n„_i(0 -f{t)) (3.35) 

where the constant matrix 91 ranges over the group SO{n — 1). 

Thus, from Theorems 13. 81 IX5l and the integrabihty properties of the Neumann system on 
rp*gn~i conclude that the phase space D C T SO{n) of the muhidimensional Veselova 
LR system with the left-invariant metric defined by (|3.15(l is almost everywhere foliated by 
(n — l)-dimensional invariant tori, on which the motion is straight-line but not uniform. 



3.4 Veselova Problem with Integrable Potentials and the Mauper- 
tuis Principle 

The Maupertuis principle. Consider a natural mechanical system on a compact Rie- 
mannian manifold {Q,ds^) with Hamiltonian h{q,p) = ^'^g^-' {q)PiPj + v{q), where g^^ is 
the inverse of the metric tensor and v{q) is a smooth potential on Q. Let By the classi- 
cal Maupertuis principle, the integral trajectories of the Hamiltonian vector field Xh with 
h{q,p) — c> maxw(g) coincide (up to a reparametrization) with the trajectories of another 
vector field Xf^j with Hamiltonian 



C^Viq) 



PiPj 



on the fixed iso-energy level £c = {h{q,p) = c} = {h'^{q,p) — 1}. Namely, on £c we have 
dh = {c — v)dh'^ (see [2])- The Hamiltonian flow of h"^ is the geodesic flow of the Jacobi 
metric dsj = (c — v{q))ds^ , which is conformally equivalent to the original metric ds^. 

The Maupertuis principle can naturally be formulated for nonholonomic systems as well. 
Suppose the distribution D is locally defined hy p = n — k independent 1-forms a* . Then 
the equations of the nonholonomic systems with Hamiltonians h and h"^ subjected to the 
constraints q G Dq are given by 

Pi ^ +^Xja^iq)i , qi = ^, i=l,...,n, (3.36) 
oqi cfpi 

Qf^J P gi^J 
= + ^/ijQ;^(g)j , q,^—, i = l,...,n. (3.37) 



1=1 



On the iso-energy level £c, the vector fields H3.36|l and H3.37|l are proportional and the 
Lagrange multipliers satisfy the relation = fii{c — v) (see |36| ) . 

One can verify that the construction goes through the Chaplygin reduction. This prop- 
erty can be used in producing non- trivial nonholonomic geodesic flows on SO{n) which, 
after the SO{n — l)-reduction, give rise to integrable systems on the sphere 5""^. 

In the case of Hamiltonian systems, under a similar reduction, the Kovalevskaya and 
Goryachev-Chaplygin integrable cases of rigid body dynamics result in integrable geodesic 
flows on S*^ that possess additional polynomial integrals of degree 4 and 3 in momenta 
respectively (see 0). 



22 



Veselova problem with potentials. Now let us go back to the n-dimensional Veselova 
problem and suppose that the n-dimensional rigid body is placed in an axisymmetric po- 
tential force field v = v{ei) (recall that {ei, . . . ,e„} are redundant coordinates on SO(n)). 
Then the equations of motion have the form 

d , . . . dv X - 

— {Ilo) = [Iw, Lo\ + - — A ei + y Xpq Bp A eg, 
at oe\ 

2<p<q<n 

ti + ioci = 0, z = 1, . . . , n, (3.38) 

togeteher with the constraints (|3.12|l . 

The potential is SO{n — l)-invariant and induces a well defined reduced potential V{q) 
on the sphere S*""^. Here V{q) = v{ei)\ei=q. The perturbed reduced system with the 
inertia tensor (|3.15() has the same Chaplygin reducing multiplier as the nonperturbed one. 
Therefore, in the new time r, the reduced system becomes a natural mechanical system on 
the sphere with the kinetic energy H3.23|l and the potential V{q). 

Let Ai < ■ ■ ■ < An- It is known, that the most general separable potentials compatible 
with the metric H3.23|) in the variables {Ai, . . . , A„_i} have the form 

where are functions of the variable Afc only (see Note that this potentials are of the 

same form as the potentials compatible with the standard metric in the same coordinates 
(e.g., see j56)'). Then, if A/j is a Laurent polynomial in the variable A^, then the potential 
II3.39|I is a Laurent polynomial in the coordinates variables qi, . . . , (see, e.g., [^ITOIIS^ ). 
In particular, the reduced Veselova problem with potential 



V{q) = ai(A-iq, q) + a2{{A-\, A-\) - {A-\, qf) + ^ 



«»+2 
1 



Ui being arbitrary constants, is completely integrable. 

Now assume that v{ei) — ai{A~^ei,ei) + a2{[A~^ei, A~^ei) — {A~^ei,eiY) and that 
the total energy is bigger than max5o(„) v. Let, as above, dsj be the left-invariant metric 
given by the inertia operator H3.15|l and introduce the Jacobi metric dsj — {c — w(ei))(isj. 
From the above considerations and the Maupertuis principle we get the following result. 

Theorem 3.9. The SO{n'-\)- reduction of the the nonholonomic geodesic flow of the metric 
ds^j with the constraints is completely integrable. The phase space y*5'"-i ig almost 

everywhere foliated by invariant (n~l) -dimensional Lagrangian tori with nonuniform quasi- 
periodic dynamics. 

The Lagrange case. In general, the operator (|3.15|) is not a physical inertia operator of 
a multidimensional rigid body. However, by taking Ai = ■ ■ ■ = /l„_i. An > Ai/2 we get 

lu^lu + ul, / = diag(/i,...,/i,/„), /i= ^ -"^'"^^ 



2detA det^ 2det^ 

In this case the system (|3.38|) represents the motion of a symmetric rigid body under the 
nonholonomic constraints. 

In the presence of the homogeneous gravitational force field in the direction ei we have 
V = Mg{C, ei), where g is the gravitational constant, M is the mass and C — (Ci, . . . , C„) 
is the position of the center of mass of the body. If the mass center is placed on the 



23 



axis of the dynamical symmetry, then v — MgCnGin and the system (|3.38|l represents a 
multidimensional version of the Lagrange top (see ^). 

In the new time t, the reduced system is completely integrable according to a non- 
commutative version of the Liouvilee theorem. Appart from the Hamiltonian function, 
there are integrals arrising from the SO{n — l)-symmetry of the system, 

QiPj " qjPi, 1 < i < j < n - 1. 

As a result, the reduced phase space T*S^~^ is foliated by two-dimensional invariant tori. 

Note that there is an another generalization of a heavy rigid body (0H|), which is based 
on the generalization of the three-dimensional Euler-Poisson equations to the Euler-Poisson 
equations on the semi-direct product so{n) x so{n). 



4 L+R Systems 

4.1 Definition and Invariant Measure of L+R Systems 

It appears that LR systems on a unimodular Lie group G can be viewed as a limit case of 
certain artificial systems on the same group, which also possess an invariant measure. The 
latter systems do not have a straightforward mechanical or geometric interpretation and 
arise as a "distortion" of a geodesic flow on G whose kinetic energy is given by a sum of a 
left- and right-invariant metrics. 



Geodesic flov^r on G with L+R metric. In addition to the nondegenerate linear opera- 
tor X defining the left-invariant metric {■,■)!, introduce a constant linear operator r° '-9^9 
defining a right-invariant metric (•, •)r on the n-dimensional compact Lie group G: for any 
vectors u, w G TgG we put {u, v)r — {ug^^ ,T°vg^^) . We take the sum of both metrics and 
consider the corresponding geodesic flow on G described by the Lagrangian 

1(^,9) = \{^,^^) + ^(.9^5"\r°gwg"^> = ^(^^,1^^) + (uj,r{g)uj), 

where T{g) = Adg-iT'^Adg and Adg is regarded as a matrix operator acting on g. 

Suppose that the total inertia operator B{g) = 2 + ^{g) is nondegenerate and positive 
definite on the whole group G. The geodesic motion on the group is described by the 
Euler-Poincare equations 

a; = [x,cj] +.g 1 — , x= — = Buj, (4.1) 

together with the kinematic equation g — g ■ oj. 

In order to find explicit expression for g^^{dl/dg), we first note that for any Y £ g, 

{Y,g-\dl/dg))=VY{l), 

where vy is the left-invariant vector field on G generated by Y. Since the metric {■,■)! is 
left-invariant, we have 

^'y{1) = ^vyUcj^Tlj)) = i(u;,radya;-Fad^ra;) = {Tcj,[Y,lu]) = (r,ad^ru;). 

As a result, g~^{dl/dg) — ad^^Fu;. 

Also, in view of the definition of F, its evolution is given by n x n matrix equation 

f = Fad^-Kad5F. (4.2) 
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Note that for compact group we have ad^ — — adt^, and f = [r,adtj]. 

Equations (|4.1(l . (|4.2|l form a cfosed system on the space g x Symm(rt) with the coordi- 
nates Wi, Tij, i < j = 1, . . . ,n. Indeed, since B is nondegenerate, the derivative to is uniquely 
defined from (|4.1|l . 

L+R systems. Now we modify equations (|4.1() by rejecting the term g^^{dl/dg). As a 
result, we obtain another system on the space g x Symm(n) 

j^{BLo)^&dlBuj, ^r = rad^+ad^r, B = l + T. (4.3) 

This is generally non a Lagrangian system, and, in contrast to equations (|4.2() . it 

possesses the "momentum" integral {Bu!,Blu). In view of the structure of the kinetic 
energy, we shall refer to the system (14 .311 as L+R system on G. 

Theorem 4.1. The L+R system j^.gl ) possesses the kinetic energy integral ^{uj,Buj) and 
an invariant measure fi dioi A • • • A aj„ A dFn A • • • A r„„ with density 



H = Vdet(J + r) . (4.4) 

Remark 4.1. The L+R systems can be also naturally considered on non-compact groups. 
Then Theorem 14 . 1 1 holds for unimodular groups as well. Recall that the group G is unimod- 
ular if tr ad^^ = 0. 

Proof of Theorem Ol First, replace -^{Buj) with Bu} + tuj. Then, using (|4.2|) and the 
identity adi^Lu — 0, we can represent equations (|4.3|l in the form 



BiJ^adlluj, f = rad„ + ad^r. (4.5) 



Using this form, we compute 

^(uj,Blj) = 2{uj,Buj) + {uj,Buj) 
dt 

= 2{Lu,8idllLu) + {uj^Tad^Lu + ad^Ftj) = 0. 

i.e., {lu,Biv) is a first integral. 

Next, divergence A of the phase flow of the system is calculated by the formula 

i<j i—1 

In view of (j4.2|l . the first sum equals J27<j + (ad(j),;i] — 0. Then we can write 



A^tiiB-'U), U,, ^ ' '\ z,j = !,...,«. 



d{a.d^Tuj)i 

^3 

As follows from the first equation in (|4.3|) , here we can put U — a,di^ + &d^I. 

In view of symmetry of B^^ , the skew symmetric part of U does not contribute to the 
expression for A. The symmetric part of U has the form 

U+ ^ \{U + U'^) = \ (ad:(6 - r) + (S - F) ad. 



25 



As a result, taking into account 14. 2|) . we obtain 

A = tr(S-i[/+) = itr (^B-^a.dlB + &d^- B-^t 
= -hviB-'t) = -ltr (B-'B). 
Now, using the unimodularity condition tr adi^ = and the well-known identity 

^detB = dctBtT{B-^B), (4.7) 



we conclude that ^ \J det B satisfies the Liouville equation ^(In/i) + A — 0, which 
establishes the theorem. 

Chaplygin's sphere. One of the best known examples of nonholonomic systems with 
an invariant measure is the celebrated Chaplygin sphere. It described a dynamically non- 
symmetric ball rolling without sliding on a horizontal plane. The center of the mass is 
assumed to be at the geometric center. Under these condition the motion is integrable 

(imini). 

It appears that a reduction of Chaplygin's sphere can be regarded as a L -I- i? system. 
Namely, the original configuration space is x 5*0(3) and the nonholonomic constraints 
define a S'£'(2)-invariant three-dimensional distribution. Then one can regard the system as 
an a Chaplygin system on the trivial bundle x 5*0(3) 5*0(3). After the R^-reduction 
we obtain a system on T5*0(3), which, written in the body frame, takes the following vector 
form 

k = Kxn, K = .m + Tna^n-ma^{n,j)-f (4.8) 
a ^ a X fl, /3 = /3xfi, = ^ x Q, 

where J, a, m, are the inertia operator, radius, and mass of the ball respectively. Next, 
is vector of the angular velocity and K is vector of the angular momentum at the contact 
point; a, /3, 7 are unit vectors forming a fixed orthonormal frame in space, 7 is assumed to be 
vertical vector. The components of these vectors can be regarded as redundant coordinates 
on 50(3). 

Equations (|4.8I) can be resolved with respect to fl to give 

9 

m = mx n + ^— (iii x 

l-ma2(7,X-i7)^ ' ^' (4.9) 

a — a X fl, $ = (3 X Q, j — j x fl. 

After usual identification of the Lie algebras {M.^, x) and (so(3), [•, •]), the system (|4.8f) 
can be seen as a L-t-R system (|4.3|l on 50(3) with left invariant metric given by the T : 
so(3) — > so(3) and right invariant degenerate operator F — ~ma?^ ® 7. According to 
Theorem l4.1l in the space (f2, 7) the above equations have an invariant measure with density 
v/det(X — ma^^ (E) 7). Up to a constant factor, it equals y^l — ma'^{'j,T ^j), the expression 
given by Chaplygin in |15| . 

Note that, in contrast to what was belived earlier, Chaplygin's sphere cannot be repre- 
sented as an LR system on the group SE{3) (see 02]). 

It is interesting that equations 1)4. 8|l are Hamiltonian with respect to a certain nonlinear 
brackets (see Borisov and Mamaev I12| V 
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4.2 The spherical Support 

The Chaplygin sphere admits an integrable generahzation on the configuration space S0{3). 
Namely, consider the motion of a dynamically nonsymmetric ball S with the unit radius 
around its fixed center. Suppose that the ball touches N arbitrary dynamically symmetric 
balls whose centers are also fixed, and there is no sliding at the contacts points. We call 
this mechanical construction the spherical support ( \22\ \2'6\ . see Figure I^TT|l . 




Figure 4.1: The spherical support 



Let g and J : — > be respectively the angular velocity vector and the inertia 
tensor of the ball 5 in a frame attached to the ball. Next, let w'^ G R^, Dk, p/c G R be the 
angular velocity, the central inertia moment and the radius of the kth peripheral ball, 7*^ 
be the fixed unit vector directed from the center C of the ball S to the point of contact 
with the kth ball, R'^ be the reaction force at this point acting on S. Then the equations of 
motion of the total mechanical system can be written in the form 

N 

Jil + nx .m = ^-f'' X R'', Dk^^ ^ -pki^ X R^, k^l,...,N, (4.10) 
fe=l 

where, as above, x denotes the standard vector product in R^^. Note that the first equation 
is taken in the moving frame, whereas the other equations are taken in a fixed frame. 

The reaction forces are due to nonholonomic constraints expressing the absence of sliding 
at the contact points. This means that velocity of the point of contact of the ball S with 
the kth ball, x 7*^ , is the same as the velocity of the corresponding point on the kth ball, 
i.e., w'^ X {—pk'l^)- Multiplying the velocities by the vectors 7*^ we obtain the constraints 
in form 

Pfe(w^ 7'=)7' - Pk^' = n-{n, 7'=)7^ k^i,...,N. 

By differentiating the constraints in the fixed frame and taking into account 7*^ = 0, 



27 



= const, we get 
and, in view of H4.10|l . 

7'= X = -:^[(f2,7^-)7'=-Si]. 
Pk 

Substituting this into the first equation in H4.1Uf) and using the fact that the time derivatives 
of r2 in the moving and the fixed frames arc the same, we obtain 

itl + Vly. JVL = -Vtl, (4.11) 

I^J-f^I, T^j:^,-^,\ (4.12) 
fe=i Pk fc^i Pk 

where I is the 3x3 identity matrix and T is the 3x3 symmetric matrix, which is fixed in 
the space. For TV > 3 and a general location of the peripheral balls, it is nondegenerate, 
hence its components can be regarded as redundant coordinates on the group S'0(3). Since 
the evolution of 7*^ in the moving frame is described by the Poisson equations t'''' = 7*^ x fj, 
from (|4.12|) we have 

r = [r,c^], (4.13) 

where lo £ so(3) is the 3x3 skew-symmetric matrix such that ujij = Eijk^k- 

Now we consider the motion of the central ball S only. As follows from H4. 1111 - 14. 13|l . 
equations of motion can be represented in the form of an L-l-R system on the group S0{3), 

k = Kxn, t=[T,uj], (4.14) 

K = {i + r)ne M^. 

Notice that from here f2 can be uniquely expressed in terms of the components of fJ, F, hence 
H4.14|) represents a closed system of differential equations. One can say that it describes the 
free rotation of a "generalized Eulcr top" , whose tensor of inertia is a sum of two components: 
one is fixed in the body and the other one is fixed in the space. 

Theorem 4.2. The spherical support system 14.14|l is integrable by the Euler-Jacobi theo- 
rem, and its generic invariant manifolds are two-dimensional tori. 

Indeed, we can put T = aa®a-\-bl3®l3-\-c^®^, where a, (3,^ are unit vectors forming a 
fixed orthonormal frame in the space and a, 6, c are some constants, which can be uniquely 
determined from ()4.12fl . Then the matrix equation in (|4.14|) can be replaced by the vector 
equations 

a = axVL, I3 = l3xn, 7 = 7x17. (4.15) 
From the form of equations (|4.14|) . H4.15|l we immediately obtain four first integrals 

iK,K), {K,a)^iIaCO,a), (K , p) ^ (1^10 , /S) , (/C, 7) = (X,c^, 7), 

Xa ^ I + 1I7 Ifc = I + 61, Tc = I + cl, 

of which any three integrals are independent. In addition, the system has trivial geometric 
integrals 

(a,a) = l, (/?,/?) = 1, (7,7) = 1, 
(a, 13)^0, (a,7)-0, (A7)=0 
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and the kinetic energy integral 

Next, according to Theorem l4.1l the system also possesses an invariant measure with density 

H = \/det(J + r) = Vdetl^l + a(a,I"^a) + b{f3,I-^f3) + c{-f,I-^j) + 

be , ^ . ac , ^ ab , ^ , abc \ 

+ d^("'^") + d^(^'^^) + d^(^'^^) + d^j ■ ^'-''^ 

This together implies the integrability by the Euler-Jacobi theorem. Notice that for the 
case of only one peripheral ball, the L+R system H4.14|l has the same form as Chaplygin's 
ball system H4.8|l . 

4.3 Limits of L+R Systems 

As mentioned above, a nonholonomic LR system on a Lie group G can be obtained as a limit 
case of a certain L+R system on this group. Indeed, suppose that the operator F : g ^ g 
defining a right-invariant metric on G is degenerate and has the form 

r = e(ai(K)a^H ha''® a''), p < n, Z) = const > 0, (4-17) 

where, as in H3.1|l . , . . . ,aP are orthonormal right-invariant vector fields a* = g^^ ■ a^ ■ g, 
= const g g, generating a right-invariant distribution D on TG. 

Now consider the L-l-R system (14.311 on the space (lu, a^, . . . , a^). In view of (|4.5|l . it can 
be represented in form 

Xtl; =X(X + r)-iad^Xw, L = Fad^ -I- ad^F. (4.18) 

Then the following theorem holds (see \2'6]). 

Theorem 4.3. 1). As e — s- oo, equations i4.18\ ) transform to the Euler-Lagrange equations 
with multipliers ^3.4)) and constraints I^S.lV . where x = Toj. 

2). The density V det Bj \fe of the invariant measure of the L+R system tends to the density 
iV.f))) of the LR system multiplied by a constant factor. 

Note that as e ^ oo, the original equations H4.3II become singular. For this reason, before 
taking the limit they must be transformed to the form (|4.18|) . 

As an illustration, consider the following L + R system on 50(3) (we use the usual vector 
notation): 

m = i{i + e-/(sj)-\mxn), 7 = 7xrj, (4.19) 

which formally coisides with the Chaplygin sphere system (|4.9|) if we set e — ma^. It can 
be easily verified that 

lim (I + e7(g)X~V)"^ =1- -7«)X~^7. 

Therefore, as e tends to infinity, the system H4.19|l transforms to the Veselova rigid body 
problem ifTTUIl . 
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